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ON GENERIC IDENTIFIABILITY OF SYMMETRIC TENSORS 
OF SUBGENERIC RANK 

LUCA CHIANTINI, GIORGIO OTTAVIANI, AND NICK VANNIEUWENHOVEN 


Abstract. We prove that the general symmetric tensor in S f!< C n '^ 1 of rank r 
is identifiable, provided that r is smaller than the generic rank. That is, its 
Waring decomposition as a sum of r powers of linear forms is unique. Only 
three exceptional cases arise, all of which were known in the literature. Our 
original contribution regards the case of cubics (d = 3), while for d > 4 we 
rely on known results on weak detectivity by Ballico, Ciliberto, Chiantini, and 
Mella. 


1. Introduction 

We denote by S' d C" +1 the space of symmetric tensors on C" +1 ; such tensors can 
be identified with homogeneous polynomials of degree d in n+1 variables, which are 
also referred to as forms. In this symmetric setting, the most natural tensor rank 
decomposition is the classical Waring decomposition, which expresses a symmetric 
tensor as a sum of powers of linear forms. Precisely, every form / £ S d C n+1 has a 
minimal expression 

r 

(1) / = £& 

i=1 

where U £ C" +1 are linear forms [25] : the minimal number of summands r is 
called the symmetric rank of /, since in the correspondence between forms and 
symmetric tensors, powers of linear forms correspond to tensors of rank 1. A natural 
question concerns the number of summands required for representing a general form 
in S d C n+1 . This problem is elementary for d = 2, which corresponds to the case 
of symmetric matrices. For d > 3, the question was answered by Alexander and 
Hirschowitz in [2]. Letting 

?n+d.\ 

( 2 ) r d>n = -^l, 

n+1 

they proved that the general / £ 5 ld C" +1 with d > 3 has rank [r^n], which is 
called the generic rank , unless the space 5 ,d C" +1 is one of the so-called defective 
cases for n = 2,3,4 and S' 3 C 5 , where the generic rank is \rd, n ] + 1- When 

the rank of a Waring decomposition is strictly smaller than r djTl , we say that this 
decomposition is of subgeneric rank. It is worth noting that, in our notation, being 
of subgeneric rank is not always equivalent to being of rank smaller than the one 
of a general tensor , because in the defective cases above a general tensor has rank 
strictly bigger than r d ,n- 
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The Alexander-Hirschowitz theorem implies that the generic tensor of subgeneric 
rank admits only a finite number of alternative Waring decompositions [251 . In 
this paper, we shall be concerned with proving that thiQ generic tensor of fixed 
subgeneric rank admits precisely one Waring decomposition, modulo permutations 
of the summands and scaling by d-roots of unity. More precisely, the main result 
of this paper is the following theorem. 

Theorem 1.1. Let d > 3. The general tensor in S d C n+1 of subgeneric rank r < 
rd,n with rd, n as In (dP has a unique Waring decomposition, i.e., it is identifiable, 
unless it is one of the following cases: 

(1) d = 6, n = 2, and r = 9; 

(2) d = 4, n = 3, and r = 8; 

(3) d = 3, n = 5, and r = 9. 

In all of these exceptional cases, there are exactly two Waring decompositions. 

The three exceptional cases were already known in the literature. The first two 
cases are classical; see Remark 4.4 in m and Remark 6.5 in m ■ The third case 
was recently found by Ranestad and Voisin; see the proof of Lemma 4.3 in m- 
A uniform treatment of these cases is presented in Proposition 12.11 Our original 
contribution establishes that there are no more exceptions to identifiability for 
cubics. 

The proof of Theorem 11.11 is based on the study of the geometric concepts of 
weak defectivity, developed in [5], and tangentially weak defectivity, developed in [B]. 
Indeed, from this point of view, the theorem can be reformulated in the following 
way, which is a result of independent interest. 

Theorem 1.2. Let d > 3, rd } n ns in m, and r < rd, n ■ Then, the common singular 
locus of the space of hypersurfaces of degree d in P" that are singular at r general 
points, consists of exactly these r points, except in the following cases: 

(1) d = 6, n = 2, and r = 9. The unique sextic plane curve singular at 9 general 
points is a double cubic, so that its singular locus is an elliptic cubic curve; 

(2) d = 4, n = 3, and r = 8. The net of quartic surfaces singular at 8 points 
consists of reducible quadrics, so that the common singular locus is the base 
locus of the pencil of quadrics through 8 general points, which is an elliptic 
normal curve of degree 4; 

(3) d = 3, n = 5, and r = 9. The common singular locus of the pencil of cubic 
4-folds singular at 9 general points is the unique elliptic normal curve of 
degree 6 through these 9 points. 

Furthermore, the above exceptional cases are the only instances where there exists 
a unique elliptic normal curve of degree n + 1 in P" through r general points. 

In this formulation, the theorem was already partially proved: the case n < 2 was 
proved by Chiantini and Ciliberto [9]; for d > 4 it was proved by Ballico 21 Theorem 
1.1]; and for d = 3 with r < rd, n — +1 it was proved by Mella m Theorem 4.1]. 

Consequently, the original contribution of this paper concerns the case of cubics, 
i.e., d = 3, which we solve completely in the subgeneric case. This answers the 
question posed in Remark 4.4 in [27j . 

We notice that Ballico [4] proved an even stronger result for d > 4. Namely, 
he showed that a general hypersurface of degree d > 4 in P” that is singular in r 

1 Since the generic tensor of fixed subgeneric rank r is an element of the r-secant variety of the 
Veronese variety Vd(P n ), which is an irreducible variety GZ], this terminology is warranted. 
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general points, is singular only at these r points (except for the exceptional cases 
(1) and (2) of Theorem 11.21) . This is equivalent to showing that the Veronese 
variety Vd(P n ) is not r-weakly defective, while our result only says that it is not 
r-tangentially weakly defective. We wonder whether the above list also gives the 
classification of all r-weakly defective Veronese varieties Ud(P"), even for d = 3. 

Symmetric tensors of general rank are not expected to admit only a finite number 
of Waring decompositions, because the expected dimension \rd,n\(n + 1) of the 
\ r d,n\ -secant variety of the Veronese variety Ud(P”) may exceed the dimension 
of the embedding space S d C n+1 . Therefore, at least a curve’s worth of alternative 
Waring decompositions of a general symmetric tensor is anticipated in these cases. 
However, if rd, n = \rd, n ~\ is integer, then a general symmetric tensor is still expected 
to admit only a finite number of Waring decompositions. The approach pursued in 
this paper, i.e., proving not tangential weak defectivity, cannot handle tensors of 
the generic rank. Other approaches need to be considered in this setting. In fact, 
Mella [28] formulated a conjecture about the cases where the expression in m is still 
expected to be unique even for general symmetric tensors. In [22] . further evidence 
for this conjecture was given; in addition, the analogous problem for nonsymmetric 
tensors was also considered. 

Even though the general symmetric tensor is not of subgeneric rank, the setting 
considered in this paper is nevertheless important in applications where one is 
mostly interested in the identifiability of symmetric tensors of subgeneric rank. 
For instance, Anandkumar, Ge, Hsu, Kakade, and Telgarsky [3| consider statistical 
parameter inference algorithms based on decomposing symmetric tensors for a wide 
class of latent variable models. The identifiability of the Waring decomposition 
then ensures that the recovered parameters, which correspond with the individual 
symmetric rank-1 terms in Waring’s decomposition, are unique, and, thus, admit an 
interpretation in the application domain. The rank of the Waring decomposition, 
in these applications, is invariably much smaller than the generic rank. As general 
sources on tensor decomposition, we refer to p~3lll6ll231l25l[34] . 

In analogy to Theorem 11.11 we mention that the results in [Q|TIJ|T2] give broad 
evidence to the analogous problem in the setting of nonsymmetric tensors, i.e., 
that a general nonsymmetric tensor of subgeneric rank admits a unique tensor rank 
decomposition, unless it is one of the exceptional cases that have already been 
proved in [HHMIDm]. 

In the proofs by induction of several theorems, we rely on the software Macaulay2 
[2T1 for proving the base cases. The two scripts we used are available as ancillary 
files in the arXiv submission of this paper. 

The content of the paper is the following. In section [2] we present a uniform 
treatment of the exceptional cases appearing in Theorems II. II and II.21 Remark l2.6l 
also discusses our initial motivation for studying the topic of this paper. Section [3] 
contains the proof of the main theorem. Thereafter, the connection between weakly 
defective varieties and the dual varieties to secant varieties, including a description 
of the dual varieties of all weakly defective examples appearing in Theorem 11.21 is 
explored in section [I] In particular, Theorem 14.21 contains the description of cubic 
hypersurfaces in P 5 that can be written as the determinant of a 3 x 3 matrix with 
linear entries. In section [5] we give an explicit criterion allowing to check if a given 
Waring decomposition is unique. This algorithm is an extension to the symmetric 
case of the one provided in 113 for general tensors. 

Acknowledgements. We want to thank C. Ciliberto, I. Domanov, J.M. Lands- 
berg, M. Mella, L. Oeding, K. Ranestad and F. Russo for useful discussions. In 
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particular, K. Ranestad pointed out the use of Gale transforms to prove Proposition 
m and informed us about Lemma 4.3 of Ell which contains a proof of the third 
exceptional case in Theorems II. II and If .21 We thank I. Domanov for pointing out 
that improved specific identifiability results can be obtained by considering reshap¬ 
ings of the tensor. The first and second author wish to thank the Simons Institute 
for the Theory of Computing in Berkeley, CA for generous support. 

2. The exceptional cases 

The following classical result shows that the values n = 2, 3, 5, which appear in 
Theorems 11.11 and ira have a special role for elliptic normal curves. 

Proposition 2.1 (Coble [T5j). A ssume that there are only finitely many elliptic 
normal curves passing through k general points in P™. Then, n = 2, 3, or 5 and, 
correspondingly, k = . In these three cases, there is a unique elliptic normal 

curve in P” passing through general points. 

Proof. Elliptic normal curves of degree (n + 1) in P" depend on (n + 1) 2 parameters, 
which is the dimension of the space of sections of the normal bundle. The passage of 
the curve through a point in P” imposes n—1 conditions, which is the codimension 
of the curve. Therefore, we may expect finitely many elliptic normal curves through 
k general points in P" only if k(n — 1) = (n + l) 2 . This implies that (n — 1) divides 
(n + l) 2 = (n — l)(n + 3) + 4, hence (n — 1) divides 4, which gives the values 
n = 2, 3, 5. Moreover, k = (n + 1 ) 2 /(n — 1). 

In case n = 2 and k = 9, the elliptic curve is a plane cubic, and it is unique. 

In case n = 3 and k = 8, an elliptic normal curve is a complete intersection of 
two quadrics. Thus, if (Qi,Q 2 ) is the pencil of quadrics through 8 general points 
Pi,... ,P8 , then C — Qi n Qi is the unique elliptic normal curve through the pf s. 

In case n = 5 and k = 9, the existence and the uniqueness of the curve was 
found by Coble [151 Theorem 19] by applying a Gale transform—see [TO] for a nice 
review—and reducing to the case n = 2 and k = 9; a modern treatment was given 
by Dolgachev [TT] Theorem 5.2]. □ 

Next, we analyze the case of cubic hypersurfaces in P 5 that are singular at 9 
general points. 

Proposition 2.2 (Veneroni [35] Section 1], Coble T5] p. 16], Room [32] Sections 
9-22], Fisher [2U] Lemma 2.9]). We have the following two results. 

(i) The 2-minors of a 3x3 matrix with linear entries on P 5 define a (sextic) 
elliptic normal curve in P 5 . 

(ii) If C is a (sextic) elliptic normal curve in P , then the variety of secant lines 
(7 2 (C) is a complete intersection of two cubic hypersurfaces on P 5 , each one 
being the determinant of a 3x3 matrix with linear entries on P 5 . 

Proof. Part (i) is well known: The curve is obtained by cutting the Segre variety 
P 2 x P 2 C P 9 , i.e., the variety of 3 x 3 matrices of rank 1, with a linear space P 5 . 
Claim (ii) follows by na Lemma 2.9]. □ 

Theorem 2.3. Let p±,... ,pg be general points in P 5 . Let C be the elliptic normal 
sextic curve through these points. A cubic that is singular at pi,...,pg contains 
a 2 (C) and is singular on C. 
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Proof. By Proposition 12.21 in the pencil of cubics containing 02 (C), the general 
element is singular along C. This pencil fills the space of cubics that are singular at 
pi,... ,pg, which is 2-dimensional by the Alexander-Hirscliowitz theorem [2]. □ 

These observations lead to a different proof of the third exceptional case in 
Theorem 11.11 

Proposition 2.4 (Ranestad-Voisin |3T[). The general tensor in S 3 C 6 of rank 9 
has exactly two Waring decompositions as sum of 9 powers of linear forms. 

Proof. In the language of [9], we have to prove that the secant order of ag (V 3 (P 5 )) 
is 2. By [S] Theorem 2.4], this is equal to the secant order of the 9-contact locus 
C, which corresponds to the third Veronese embedding of an elliptic normal sextic 
curve in IP 5 , by Theorem l2.3l Thus, C is an elliptic curve of degree 18 in P 17 , whose 
secant order is 2 by 0 Proposition 5.2]. □ 

Corollary 2.5. Let n = 2, 3, or 5. The unique elliptic normal curve that is 
mentioned in Proposition \2.1l which passes through general points in P n , 

can he constructed as the singular locus of a general hypersurface of degree 
that is singular in the points. 

Proof. The cases n — 2, 3 have already been considered in the proof of Proposition 
cm The case n = 5 follows from Theorem 12. 31 □ 

Remark 2.6. It was only at the completion of this manuscript that we became 
aware of Ranestad and Voisin’s proof of the third exceptional case that appears in 
Theorems 0 and O Our initial motivation for studying this problem arose be¬ 
cause the third case was unexpectedly — in our minds—suggested by a computational 
analysis performed by the third author, who ran the algorithm that we present in 
section [5] for hypersurfaces of degree d in P n singular at the maximal number of 
random points, i.e., r = rd, n — 1 with rd, n as in 0 ), for all reasonably small values 
of d, n. It took a while to realize what happened, because this third case was missing 
in V261 Theorem 6.1.2]. Actually, Theorem. 6.1.2 of f26] only intended to collect pre¬ 
vious results by Ballico m Ciliberto and Chiantini m, and Mella \27^28f . which are 
individually correct. The second author takes the responsibility to have first over¬ 
looked the assumption d > 4 in summarizing and reporting the results of 0 |H7f . For 
d = 3, Theorem \l.l\ was known with the additional assumption r < r^ yn — ^—I- 1; 
see \27\ Theorem ].l]. From the theoretical proof that we present in section\3 ] we 
can conclude that the third case was the last exception. Therefore, Theorem 6.1.2 
in \26] remains true if the third case ( k , d, n) = (9, 5,3) is added to the list of excep¬ 
tions. Exactly the same remark applies to the formulation of Theorem 2.3 in fW] 
and Theorem 12.3.4-3 in We informed the coauthors of !26i29f of the problem, 
and they accepted the above conclusion. 


3. Cubics singular at the maximum number of points. 


We turn our attention to the proof of Theorem 11.11 in the case of cubics, i.e., 
d = 3. Given n , we define 


r ( n s 3 ) 1 


~{n + 3)(ro + 2)' 

n + 1 


6 


it is the generic rank for cubic polynomials for n / 4. In other words, a cubic 
polynomial on P n singular at k n general points vanishes identically for n 7 I 4 
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[2j. Some elementary algebra shows that k n = ( n+ 3 H ra + 2 ) if n ^ 2 mod 3, while 
k n = (n+3Kn+2) + 2 = (n+4Kn+l) + X jf „ = 2 mo d 3. 

For the sake of future reference, let us state explicitly the following consequence 
of the Alexander-Hirschowitz theorem [2] . 

Theorem 3.1 (Alexander-Hirschowitz [5j). The space of cubic hypersurfaces on 
P" that are singular at k n — 1 general points has dimension 

(i) n+1 if n ^ 2 mod 3, or 

(ii) if n = 2 mod 3. 

In addition, the space of cubic hypersurfaces on P n , n / 4, that are singular at k n 
general points is empty. 

To complete the proof of the Theorems 11.11 and 11.21 it remains to show the 
following result, which refines Theorem 13.II 

Theorem 3.2. The space of cubic hypersurfaces in P" that are singular at k n — 1 
general points has dimension 

(i) n+1 if n ^ 2 mod 3, or 

(ii) if n = 2 mod 3, 

and, in addition, its common singular locus consists only of these k n — 1 points, 
provided that n ^ 5. 

In order to prove Theorem 13. 21 we may assume n > 6 , since the cases with n < 4 
(as well as the case of cubics in P 5 singular at 8 points) can be checked separately 
using the approach described in section [5] The outline of our proof of Theorem 
13.21 is as follows. In section IQ we will prove case (i) by induction on subspaces 
of codimension 3, adopting an approach that is mainly inspired by [7, Section 
5], where an alternative proof of Theorem 13.11 was given. To prove case (ii), the 
aforementioned technique needs a modification. We will construct an inductive 
proof on subspaces of codimension 3 and 4; in the inductive step, we will rely, 
additionally, on the argument of case (i). This strategy will be presented in section 

In the rest of this section, if 5 is a set of simple points in P" and P C P" is 
a linear subspace, we denote with Is,p(d) the space of degree d polynomials in P 
vanishing at all of the points in S. Moreover, if X is a a set of double (singular) 
points, we denote by /xus,p(d) the space of degree d polynomials in P vanishing 
on all of the points in iSUX and whose derivatives vanish on all of the points in X. 
The notation L = (xi... Xi+d) denotes the subspace of codimension d + 1 whose 
ideal is (xi,x i+1 ,..., x i+d ). 

3.1. Proof of Theorem 13.21 (i) by induction on codimension 3. We start by 
proving three auxiliary results. 

Proposition 3.3. Let n > 6 , and let L , M, N C P" be general subspaces of codi¬ 
mension 3. Let U, respectively mi, with i = 1,2,3 be three general points on L, 
respectively M. Let ni with i = 1,2 be two general points on N. Then, the space 
of cubic hypersurfaces in P" that contain L U M U N and that are singular at the 
eight points X = {Zi, I 2 , 13 , m±, m 2 , m 3 , ni, 712 } has dimension 3. Furthermore, the 
common singular locus is contained in LU M U N. 

Proof. The base cases n = 6 , 7 and 8 can be proved with the Macaulay2 script 
generic-identifiability.m2 that is provided as an ancillary file to the arXiv 
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version of this paper. Using this software, we may compute the following dimen¬ 
sions: 

dim/xuLuJVuM,pe(2) = 0, dim Ixuluivum,?** (3) = 3, and 

dim/ X uLuivuM,P 7 (2) = 0, dim/ X uLuivuM,P 7 (3) = 3, 

so that the claim on the codimension follows. The code also proves the statement 
about the singular locus. 

For n > 9, the statement follows by induction on n. Indeed, we may choose 
coordinates such that L = ( X 0 ...X 2 ), M = (X 3 ...X 5 ), N = ( xe...xg ). In this 
setting it is clear that there are no quadrics that contain IUMU1V, and moreover 
every cubic containing L U M U IV is a cone with vertex in L fl M fl IV. Thus, for 
a general hyperplane H C P ra , the Castelnuovo sequence (see (7] Equation (1)]) 
induces an inclusion 

0 - >lLUMUN,P n (3) - >I(LUMUN)r\H,H(3)- 

Hence, if we specialize the eight points to the hyperplane H , we get an inclusion 
0 -t/xuLUMUIV,P"(3)->-d(xuLUMUJV)n7f,7f(3)- 

Then, our statement follows by induction. The singular locus is a cone with vertex 
L n M fl N over the singular locus of the base case n = 8. □ 

Remark 3.4. Following the output of the software for the case n = 8 , we can 
guess the common singular locus of the cubic hypersurfaces in Ixulumuiv,p ,i (3), 
for n > 8. It turns out that in some examples—but we believe in general—it is 
given by the union of the three linear subspaces L n M, L fl IV, M fl N and by 8 
linear subspaces of codimension 7, each containing one of the 8 points, and three of 
them contained in L, three of them contained in M, and two of them contained in 
N. 

Proposition 3.5. Let n > 5, and let L , M C P" be subspaces of codimension 
three. Let li, respectively mi, with i = 1,... ,n — 2 be general points on L, re¬ 
spectively M. Let pi,P 2 £ P" be general points. Then, the space of cubic hy¬ 
persurfaces in P” containing L U M and singular along the set of 2n — 2 points 
X = {7i, I 2 , ■ ■ ■, l n - 2 , mi, m 2 , ■ ■ ■, m n - 2 ,Pi,P 2 } has dimension n + 1. Its common 
singular locus contains the linear space L fl M and is 0-dimensional at the points 
Pi and P 2 ■ 

Proof. The base cases n = 5, 6, and 7 can be proved with the Macaulay2 script 
generic-identifiability.m2. Running the software, we find the following di¬ 
mensions 

dim I X uiuM,P 6 (2) = 0, dim/ X ulua7,p 5 (3) = 6, 

dim/ X uLuAf,P 6 (2) = 0, dim/ X ulua7,p 6 (3) = 7, and 

dimIxuiuAf,P 7 (2) = 0, dim/ X ulua7,p 7 (3) = 8. 

These values indeed correspond to the claimed dimensions. 

For n > 8, the statement follows by induction from n — 3 to n. Indeed, given a 
third general subspace N of codimension 3, we get the exact sequence 

0 >lL\JMUN,P n (3) >lLUM,P n (3) >I(LUM)r\N,N{ 3), 

where the dimensions of the three spaces in the sequence are respectively 27, 9 [n — 
1), and 9(n — 4). Let us specialize n — 5 of the points k £ L to L fl N, n — 5 of 
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the points rrii € M to M n N, and the two points pi,P 2 to N. Then, we obtain a 
sequence 

0 -t JxuiUMUiV.P’ 1 (3) TxuLUMjP" (3) t/(xu£UM)rW,iv(3), 

where the trace (XULU M) n N satisfies the assumptions on iV = P"~ 3 , so that 
we can apply the induction. Notice that the residual (left) space satisfies the hy¬ 
potheses of Proposition 13.31 and has dimension 3. Since the common singular locus 
of the cubics containing LLlM and singular at X must be contained in the common 
singular locus of the leftmost 3-dimensional space, it follows by Proposition l3.3l that 
its components through p\ and P 2 must be contained in N. After the degeneration, 
the space of cubics /xuium.P” (3) still has dimension at most 3 + (n — 2) = n + lby 
induction. Hence, by semicontinuity it follows that its dimension is indeed equal to 
7 i+l. The common singular locus cannot be positive dimensional at points p\ and 
P 2 , because otherwise it should be of positive dimension in the trace (right space), 
where by induction we know that it is O-dimensional. □ 

Proposition 3.6. Let n > 6 , and let L C P" be a subspace of codimension 3. If 
n ^ 2 mod 3, then the space of cubic hypersurfaces in P n that contain L and that 
are singular at general points I, £ f and at n general points pi € P" has 

dimension n + 1. Moreover, its common singular locus is O-dimensional at the n 
points Pi. 

Proof. The statement can be checked for n = 6 , 7 with the Macaulay2 script 
generic-identifiability.m 2 . 

Let n > 9 and n jk 2 mod 3. Consider the sequence 

0—t/LuM,P"(3)— >Il, P n (3)— >lLr\M,M(3), 

where M is a general subspace of codimension 3. Denoting by X the union of the 
double points supported at the points U and pi, we get 

0-dxuLUM,P n (3)-iTxuL,P"(3)-t/(xui)nM,Af (3)- 

We specialize ( w ~ 3 H rt ~ 4 ) 0 f th e points U to L D M and n — 2 of the points pi to 
M. We can assume that at least one point pi that is a contained in a positive 
dimensional component of the singular locus is not specialized. Thus, we left n — 2 
general points on L and 2 general points in P". Let us note that we cannot apply 
induction from n — 3 to n to determine the dimension of the right space—contrary 
to the strategy that was employed in the proof of the foregoing propositions in this 
section -because then we would have to specialize n —3 of the points pi to M, hereby 
losing control over the singular locus. Instead, we note that we can immediately 
use Proposition 5.4 of [7j (in P”~ 3 ) on the trace (right space); it turns out to be 
empty. On the residual (left space), Proposition 13.51 can be invoked, proving that it 
has dimension n + 1. Now, if the singular locus would have a positive dimensional 
component, then, since the dimension of the space of cubics is constant along the 
specialization (it equals n + 1), we would get a deformation of the singular locus, 
which should be of positive dimension at every point. This, however, contradicts 
Proposition [331 hereby concluding the proof. □ 

We are now ready to prove the first part of Theorem [33] 

Proof of Theorem \3.2[ case (i). We fix a linear subspace L C P" of codimension 3 
and consider the exact sequence 
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where Sp»(3) is the space of cubic polynomials on P" and the quotient space Sl{ 3) 
is isomorphic to the space of cubic polynomials on L. Then, we specialize to L as 
many points as possible in such a way that the trace with respect to L imposes 
independent conditions on the cubics of L. To be precise, we have k n — 1 = 
(n+3)(n+2) _ ^ (j ou bi e p 0 j n ts and we specialize k n -3 = nflLzl1 G f them to L, leaving 
n points outside. Then, the result follows from Theorem 5.1 of [7] on the trace 
(right space), which turns out to be empty, and by Proposition 13.61 on the residual 
(left space), which has dimension n + 1. If the contact locus has positive dimension, 
then, since the dimension of the space of cubics is constant and equal to (n + 
1 ) in the degeneration, we would get a deformation of the singular locus with a 
positive dimension at every point, contradicting Proposition 13.61 and concluding 
the proof. □ 

3.2. Proof of Theorem 13.21 (ii) by induction on codimension 3 and 4. 

For proving the second case in Theorem 13.21 we need to introduce several other 
auxiliary results on configurations that involve subspaces of codimension three and 
four. These configurations are covered in Propositions 13.71 through 13.121 

3.2.1. Codimension 4, 3, 3. 

Proposition 3.7. Let n > 6, and let L, M, N C P" be general subspaces of codi¬ 
mension 4, 3, and 3, respectively. Let h,l 2 ,h be general points on L. Let mi, 
respectively, ni with i = 1,..., 4 be four general points on M , respectively N . Then, 
the space of cubic hypersurfaces in P n that contain LU M U N and are singular at 
the 11 points X = {Zi, I 2 , 13 , mi, ..., 7714, n\, ..., 714} is empty. 

Proof. The proof is similar to the proof of Proposition 13.31 The Macaulay2 code 
proves the base cases n = 6, 7, 8 and 9. For n > 9, we may choose coordinates such 
that L = {xq ... X 3 ), M = (X 4 ... xg), and N = (ay ... xg). Then, the statement 
follows by induction on n. Indeed, as in the proof of Proposition 13.31 the space 
^LuMuAr,P n (2) is empty, thus for a general hyperplane H C P”, the Castelnuovo 
sequence induces an embedding 

0 >IluMUN,P™(3) - >I(LUMUN)nH,H(3), 

and, moreover, every cubic in the left space is a cone with vertex at L fl M l~l N. 
Hence, by specializing the 11 points to the hyperplane H, we get: 

0 >IxuLUMUIV,P"(3)- >I(X.ULUMUN)r\H,H{3)- 

Then, the statement follows by induction. □ 

Proposition 3.8. Let n > 7, let n = 1 mod 3, and let L,M C P" be subspaces of 
codimension 4 and 3, respectively. Let U with i = 1,..., n — 3 be general points on 
L. Let mt with 7 = 1,..., 4 ”~ ln be general points on M. Then, the space of cubic 
hypersurfaces in P n that contain Lull and are singular at all the L’s and mi’s 
and at four general points pi,P 2 1 PSiPi & P n , is empty. 

Proof. The Macaulay2 script proves the base case n = 7. 

For n = 3k + 1 with k > 3, the statement follows by induction from n — 3 to 
n. Indeed, given a third general subspace N of codimension 3, we get the exact 
sequence 

0 —t Ilumun ,p n (3)—>Ilum,p» (3)— ^I(lum)c\n,n (3), 

where the dimensions of the three spaces in the sequence are respectively 36, 12n—18 
and 12n — 54. Let X denote the union of the double points supported at the pf s, 
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L’s and mi’s. Assume that we specialize n — 6 of the points k G L to L n TV, 4 ”~ 22 
of the points rrii € M to M f~l TV, and the four points p\,... ,p 4 to TV. Then, we 
obtain a sequence 

0-tdxuLUMUiV,P™(3)-dxuLUM,P" (3)->-f(XULUM)niV,iv(3), 

where the trace (XULUM) nTV satisfies the assumptions on TV = P™ -3 , so that we 
can apply induction. Then, we may conclude, as the residual (left space) satisfies 
the hypotheses of Proposition 13.71 and, consequently, it is empty. □ 

Proposition 3.9. Let n > 7, n = 1 mod 3, and L C P n be a subspace of codi¬ 
mension four. Then, the space of cubic hypersurfaces in P n that are singular at 
k n -4 = ( n ~ 1 h n ~ 2 l general points h on L (and, thus, contain L, by Theorem \3.1\) 
and at general points pi G P n is empty. 

Proof. The Macaulay2 script generic-identif iability .m2 proves the case n = 7. 
For n = 3k + 1 with k > 3, the statement follows by the sequence 

0—>-/lum,P"(3)— >Il, P n (3)—3), 

where M is a general subspace of codimension 3. If we denote by X the union of 
the double points supported at the points U and pi, then we get the sequence 

0-t/xuLUM.P" (3)-t/xuL,P"(3)-H(xuL)nM,M(3). 

Then, we specialize k n -7 of the points k to L fl M and 4 ”~ 10 of the points Pi 
to M. The trace (right space) contains exactly fc n _3 double points and turns out 
to be empty by induction. Thus, there remain n — 3 general points on L and 4 
general points on P”; we can then use Proposition 13.81 on the residual (left space) 
to conclude. □ 

3.2.2. Codimension f, f, 3. 

Proposition 3.10. Let n > 8 , and let L,M,N C P" be general subspaces of 
codimension respectively 4, 4, and 3. Let li, respectively nii, with i = 1,... ,4 be 
general points on L, respectively M. Finally, let ni with i = 1,...,5 be general 
points on N. Then, the space of cubic hypersurfaces in P" that contain L U M U N 
and that are singular at the 13 points X = {Zi,..., Z 4 , mi,..., 7714 , m,..., 715 } has 
dimension 1. In other words, there is a unique cubic hypersurface W through L U 
MU JV and singular at X. Furthermore, the singular locus of W is contained in 
LUMUN. 

Proof. The proof is similar to the proof of Proposition 13.31 The Macaulay2 code 
generic-identif iability .m 2 proves the base cases n = 8 , 9, and 10. 

For n > 11, we may choose coordinates such that L = {xq ... xf), M = (24 ... X 7 ), 
and N = (Xg ■ •. X 10 ); then, the statement follows by induction on n. Indeed, as 
in the proofs of Proposition 13.31 and Proposition 13.71 we let H C P" be a general 
hyperplane, so that the Castelnuovo sequence induces the inclusion 

0 -diuMUiV l P n (3)- >I(LUMUN)r\H,H( 3), 

because the space Ilumun.p "(2) is empty. Hence, by specializing the 13 points on 
the hyperplane H , we get an exact sequence: 

0 ^-ZxuLUMUiVjP 71 (3) >I(JQULUM\JN)nH,H (3). 

Now the statement follows by induction. □ 
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Proposition 3.11. Let n > 8, n = 2 mod 3, and L,M C P" be subspaces of 
codimension 4. Let U and mi, where i = 1,, 4 ”" 14 , be general points on L and 
M, respectively. Then, the space of cubic hypersurfaces in P” that contain III M 
and are singular at the 8 ”~ 28 points h,mi, i = 1,..., 4 "~ 14 , and at an additional 
set of five general points pi G P”, i = 1,..., 5, has dimension Furthermore, 
its common singular locus, which contains the linear space LC\M, is 0 -dimensional 
at each of the points p \,..., ps. 

Proof. The case n = 8 is handled in the generic-identif iability .m2 Macaulay2 
script. 

For n = 3k + 2 with k > 3, the statement follows by induction on k. Given a 
third general subspace N of codimension 3, we get the exact sequence 

0—diuMuJV,p»(3)—tJz,uM,P n (3)—>-f(LuM)niv,iv(3), 

where the dimensions of the three spaces in the sequence are respectively 48, 16(n — 
2) and 16(n — 5). 

Let X denote the union of the double points supported at p\,... ,p$, h and mi 
with i = 1,..., 4ra ~ 14 . Then, we specialize 4n f 26 of the points Zj G L to L (~l N, 
4ra ~ 26 of the points mi G M to M fl N, and the points pi,... ,p 5 to N. We thus 
obtain a sequence 

0—^xulumujv,P"(3)—>Txulum,p™( 3)—»-f(xuiuM)niv,iv(3), 

where the trace (XUlU M) n N satisfies the assumptions on TV = P n ~ 3 , so that 
we can apply induction. Then, the residual (left space) satisfies the hypotheses 
of Proposition 13.101 and has dimension one. Moreover, the common singular locus 
has to be contained in the common singular locus of the left 1-dimensional space. 
After the degeneration, the space /xulum,P"( 3) still has dimension less than or 
equal to 1 + by induction, and, therefore, its dimension equals by 

semicontinuity. The common singular locus cannot be positive dimensional at the 
points p\,... ,P 5 , because otherwise it should be positive dimensional in the trace 
(right space), while we know that it is 0-dimensional there by induction. □ 

Proposition 3.12. Let n > 8, n = 2 mod 3, and L C P" be a subspace of 
codimension 4. Then, the space of cubic hypersurfaces in P" that contain L and 
are singular at k n -4 = ( n ~ 1 h Ti ~ 2 ) general points U G L and at 4 ” 4 ~ 1 general points 
Pi G P™ has dimension . Furthermore, its singular locus is of dimension 0 at 
all of the points pi. 

Proof. The statement follows by the sequence 

0—diuM.P" (3) — >Il ,P" (3) —(3), 

where M is a general subspace of codimension 4. Denoting by X the union of the 
double points supported at the points If s and pf s, we get 

0 dxuLUM,P n (3) t/xuL,P"(3) bf(xui)nM,M(3). 

Suppose that the singular locus would be of positive dimension at one of the 
pf s, say, at q = pj. Then, we specialize of the points k to L fl M and 4 ”~ 14 
of the points {pi}i^ij to M. The trace (right space) contains exactly k n -4 double 
points and is empty because of Proposition [T9] There remain 4 ”~ 14 general points 
on L and 5 general points—one of which is q —on P n . We can use Proposition 
13.111 on the residual (left space), which has dimension By assumption, the 

singular locus has positive dimension at q; however, the dimension of the space of 
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cubics is constant and equal to through the degeneration, so that we have a 
deformation of the singular locus, which is of positive dimension at every point. In 
particular, the singular locus will be of positive dimension at all the points that 
were not specialized to M, hereby contradicting Proposition 13.111 We conclude 
that our initial assumption must have been false, so that no general points p 3 can 
exist where the singular locus is of positive dimension. □ 

Proof of Theorem \3.2l part (ii). We fix a codimension four linear subspace L CP" 
and we use the exact sequence 

0 — >Il, p "( 3 )— >Spn(3 )— >Sl( 3 ), 

where, as above, <Spn(3) is the space of cubic polynomials on P” and the quotient 
space Sl{ 3) is isomorphic to the space of cubic polynomials on L. We specialize 
k n -4 points on L , leaving 4 ” 3 f 1 points outside. Then, the result follows from The¬ 
orem 5.1 of S.7j on the trace (right space), which turns out to be empty and by 
Proposition 13.121 on the residual (left space). If the contact locus would have a 
positive dimension, then, since the dimension of the space of cubics is constant and 
equal to in the degeneration, we would get a deformation of the singular locus, 
which should be of positive dimension at every point; however, this contradicts 
Proposition 13. 121 hereby concluding the proof. □ 

4. Dual varieties to the relevant secant varieties 

Denote by T x X the tangent space to the projective variety X C P” at the point 
x £ X. Following the notation of [9] we say that X is not k-weakly defective if the 
general hyperplane H containing the linear span of the tangent spaces at k general 
points x\,.. ., Xk £ X, i.e., (T Xl X, ..., T Xk X) C H , is tangent to X only at finitely 
many points. This is equivalent with saying that the fc-contact locus with respect 
to Xi,... ,Xk and H is zero-dimensional. 

For any projective variety X, we will denote by X v the dual variety to X. Note 
that the dual of the secant variety ijfc(i; [ j(P n )) v contains the points correspond¬ 
ing to hypersurfaces of degree d in P” with k general singular points, and it has 
codimension > k, where k is the expected value for the codimension. 

Proposition 4.1. Let X C P N and let <Jk{X) be the k-secant variety of X. Then, 
the following are equivalent: 

(i) the general hyperplane H containing (T Xl X, ..., X Xk X) for general aq,..., Xk 
is tangent to X only at aq ,...,Xk, i.e., the k-contact locus with respect to 
x\ ,... ,Xk and H consists exactly of the points aq, ... , Xk, 

(ii) X is not k-weakly defective, and 

(Hi) dim [crfc(A')] = N—k, that is a general hyperplane tangent to <7 k{X) is tangent 
along a linear space of projective dimension k — 1 . 

Proof, (i) -£=>• (ii) follows from [ 8 ] Theorem 1.4], (ii) (iii) follows from Ter- 
racini’s Lemma. □ 

It is interesting to describe the dual varieties of crfc(ud(P")) in the exceptional 
cases of Theorems D and II.21 They have dimension smaller than expected. 

Theorem 4.2. The following dual varieties correspond to the exceptional cases 
appearing in Theorems \l.l\ and \1.2\ 

(i) o 9 (ug (P^)) v contains the plane sextics which are double cubics. It has codi¬ 
mension 18. 









IDENTIFIABILITY OF SYMMETRIC TENSORS 


13 


(ii) ct 8 (w 4 (IP 3 )) v contains the quartic surfaces which are reducible in a pair of 
quadrics. It has codimension 16. 

(Hi) ag(v 3 (P' 5 )) v contains the cubic 4- folds which can be written as the determinant 
of a 3x3 matrix with linear entries. It has codimension 18. 

To compute the dimension in third case, note that the Hilbert scheme of elliptic 
normal sextic curves in P 5 has dimension 36. So the cubic hypersurfaces coming 
from this construction have dimension 37, and 37 + 18 = 55 = (®) — 1. 

We remark that the defective Veronese varieties according to the classification 
of Alexander and Hirschowitz [2] (see [30] for the equations of the defective secant 
varieties) yield the following dual varieties 

(i) u" n ( n _(_ 3 )/ 2 {va (P n )) v , for n = 2,3,4, contains quartic hypersurfaces which are 
double quadrics. It has codimension ( n ^ 2 ) n -^ J -- 

(ii) 07 (i; 3 (P 4 )) v contains cubic 3-folds which can be written as the determinant of 
a 3 x 3 symmetric matrix with linear entries. It has codimension 13. Indeed, 
it is birational to the Hilbert scheme of quartic rational normal curves which 
has dimension 21. 

5. Specific identifiability of symmetric tensors 

While the generic symmetric tensor of subgeneric rank is expected to admit a 
unique Waring decomposition, specific tensors, whose Waring decomposition is as¬ 
sumed to be known, may admit multiple decompositions. We proceed by presenting 
an approach for certifying specific identifiability of symmetric tensors of small rank 
by checking not tangential weak detectivity of the r-secant variety of a Veronese 
variety in the given point. The strategy is an adaption of the algorithm from m 
to the setting of identifiability with respect to the Veronese variety V = Vd(P"). 
As such, the presented condition will only be a sufficient condition; that is, if the 
criterion does not apply, then the outcome of the test is inconclusive. On the other 
hand, if the criterion applies, then the given input tensor is r-identifiable and of 
symmetric rank r. Throughout this section, it is assumed that we are handed a 
Waring decomposition 

p = pi H- +p r £ oy(V) C S d C n+1 , 

wherein the point pi = a.f d £ V is the degree d Veronese embedding of the vector 
a i £ C n+1 . In other words, we know the points pi appearing in the decomposition. 
The goal only consists of certifying that p is r-identifiable—the decomposition(s) 
are not sought. To this end, the strategy in m suggests a two-step procedure: 
Prove that p is a smooth point, and verify the Hessian criterion. In principle, the 
method can be applied for all tensors of subgeneric rank, however, in practice the 
range of applicability of the algorithm in section 15.31 is restricted by the lack of 
good techniques for certifying smoothness. For this reason, the results of [18] may 
apply in a wider range than the results we present, by combining reshapings of 
higher-order tensors into tensors of order three with Proposition 1.32 in [Tj3], as we 
were kindly informed by its authors. We will nevertheless present an example of 
an identifiable tensor in S' 3 C 7 of rank 10 whose identifiability cannot be proved by 
the state-of-the-art specific identifiability criteria from the literature. 

It is important to stress that we discuss the general setting of degree d > 3 
Veronese embeddings. We will restrict our attention to nondefective r-secants of 
V, because identifiability will not hold for general tensors on a defective r-secant 
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variety. This is the interesting setting, because the Alexander- Hirschowitz theorem 
[ 2 ] stipulates that most ay (V) are nondefective. 

5.1. The Hessian criterion. We recall the main proposition from mi and adapt 
it to the present context of symmetric tensors. 

Lemma 5.1 (Sufficient condition for specific identifiability). Let V = u<j(P™) be a 
nondefective Veronese variety, and let r < |— 1 with rd, n as in fH). Assume 
that we are given a nonsingular point 

P = Pi+P 2 ~\ - \~Pr G oy(V). 

If the linear span of the tangent spaces to V at the pi’s, i.e., 

M = (T Pl V,..., T Pr V), 

has the expected dimension, i.e., r(n + 1 ), and if, in addition, the r-tangential 
contact locus 

C r = {p G V | T p V C M} C V 

is zero-dimensional at every p\,p 2 , ■ ■ ■ ,p r , thenp is r-identifiable, r is its symmetric 
rank, and p = ^i_i Pi is its unique decomposition. 

Proof. The proof is obtained by repeating the proofs of m Lemma 4.3, Lemma 
4.4, and Theorem 4.5], therein substituting the Segre variety with the Veronese 
variety V. We present a simplification of the proof of mi Theorem 4.5]. There is 
an open neighborhood of p = p± + • • • +p r consisting of points for which smoothness, 
Terracini’s lemma [25]; and the absence of a contact locus will hold. The variety 
must thus be generally identifiable, hence the projection 7r onto the first factor 
of the usual abstract secant variety Aa r (V) is a birational morphism. After we 
find another decomposition p = we get that the fiber 7 T~ 1 (p) contains 

the two points (p, (p i,.. . ,p r )) and (p, (q±, ..., q r ))- Terracini’s lemma implies that 
the connected component of the fiber passing through (p, (pi,... ,p r )) cannot be 
positive dimensional, hence, it contains just this unique point. Since 7r is birational 
and p is a smooth point of oy(V) = ir(Aa r (V)), we have a contradiction with 
Zariski’s Main Theorem. □ 

Remark 5.2. We note a minor omission in the formulation of Theorem 4-5 in m > 
where we forgot to include the condition that M should be of the expected dimension. 
It is clear from the proof of aforementioned theorem that this condition must hold, 
as can be understood from the invocation of JT3. Lemma 4-3]. 

Remark 5.3. If one chooses r random points pi £ v ( [(W n ), then Lemma \5.1\ may 
be invoked to prove generic r-identifiability. In this way, one can handle the cases 
773 (P 2 ), V 3 (P 3 ) and i> 3 (P 4 ), which were not covered by the proof in section [H The 
case ^(P 1 ) is trivial, because there is only one point, which is naturally identifiable. 

For practically verifying Lemma 15.11 we need a sufficiently explicit description 
of the r-contact locus. This is obtained as follows. Interpreting a point ft e V as a 
power of a linear form, say 

Pi — (ao,;£o + T * * ■ T , 

where {a7i }"_ 0 is a basis of P", it follows immediately that the tangent space is 
given by 

T Pi V — (xq(oq^xq-{- ai^Xi -{-• * *-f ~a n ^iX n ') ,..., x n (ao^iXo-\-ai^iXi-\-‘ * *-f ~a n ^x n ') ). 
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If we choose the standard monomial basis {xi ± Xi 2 ■ ■ ■ Xi d }o<i 1 <i 2 <--<i d <n f° r Ud(P"), 
then this tangent space can be represented in a straightforward manner as a x 
(n + 1) matrix of constants, say T], The Cartesian equations of M may then be 
constructed by computing the kernel of the matrix T = [Ti t 2 ■■■ tv] . The number 
of such equations should be precisely £ = ( n ~j i d ) — r(n + 1); otherwise, the first 
condition in Lemma O concerning the dimension of the tangent space would be 
violated. Let us denote the Cartesian equations as 

n n n 

(5) Ql (x) ^ ^ ( * * * ^ ( ^(ii ,i2,.. ’ “Cl *^22 * * ' ^ 1, 2, . . . , 

il=0i 2 =il id=id-1 

where the vector k; = [fc(i ll i 2l ... I i (J ).z]o<i 1 <i 2 < --<id<»T. is the Zth basis vector of the ker¬ 
nel of the matrix of constants T. For imposing that a point p = (a oxo + ■ ■ ■ + a n x n ) d 
is contained in M, it should obey the Cartesian equations, i.e., qi(ao, <zi,..., a n ) = 0 
for all / = 1,2,..., i. That is, the equations |3]) define the ideal-theoretic equations 
for V fl M. It similarly follows that deriving the equations m with respect to 
Xq, xi, ..., x n and substituting xo, Xi, ..., x n for, respectively, a o, a±,... ,a n results 
in the ideal-theoretic equations of the intersection C r = T p V fl M; naturally, the 
ao, cii,. •., a n should be treated as new variables. The number of equations thus 
constructed equals £{n + 1). To determine that C r is zero-dimensional at each pi, it 
suffices to verify that the codimension of the tangent space, i.e., the derivative of 
the equations of the ideal, is n at each of the pi s. This tangent space can be repre¬ 
sented by a matrix H of size (n-l- 1) x £(n+ 1), which contains only constants when 
it is evaluated at one of the pi ’s. The rank of H coincides with the dimension of the 
contact locus and can be computed using simple linear algebra. As was remarked 
in m, H can be interpreted as a “stacked Hessian” matrix H = [h 1 h 2 ... h 1 }, 
wherein H k is the Hessian matrix of partial derivatives 

d 2 I" 

qk (X 0 , Xl, . . . , Xn) , 

OXjOXi Jij=0 

this is the reason why we call the above approach of verifying Lemma 15.11 the 
Hessian criterion. 

A computer implementation of the Hessian criterion in Macaulay2 is included 
in the specif ic-identif iability .m2 file that accompanies the arXiv version of 
this article. 

5.2. The smoothness criterion. The Hessian criterion in Lemma lSTl mav only be 
applied to smooth points of ay(V). One approach for proving smoothness consists 
of verifying that the local equations of ay(V) are of the expected degree. Such 
equations are known in the case when the number of terms r in the symmetric 
decomposition is sufficiently small. A standard nontrivial set of local equations is 
generated by the (r+l)-minors of the usual symmetric flattenings; see [231 Theorem 
7.3.3.3] and [23j Theorems 4.5A and 4.10A]. For Veronese embeddings of odd 
degree, the Young flattenings from [251 Section 4] apply in a wider range than the 
standard symmetric flattenings; however, they are more involved to explain and 
implement. Our discussion will focus on the simple symmetric flattenings, which 
can still handle a respectable number of cases for Veronese embeddings of degree at 
least four. For degree three Veronese embeddings, the Young flattenings that were 
described in [25] should be employed. 

The strategy that we present for proving that p corresponds to a smooth point 
is, essentially, based on [25] Theorems 4.10A and 4.5A] and [23] Theorem 7.3.3.3] 


H 
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and consists of obtaining local equations of the r-secant variety ay (V) in p. Crucial 
to this approach are the symmetric flattenings , which, we recall, may be defined as 
follows. Let p E S d C n+1 , then we can define the map 

cj? k : S k (C n+1 )* -> s d ~ k C n+1 

d k 

1 ^ "a a a ' 

OXi-y C/Xi 2 * ’ * ox ik 

We have the following. 

Lemma 5.4 (Sufficient condition for smoothness). Let V = i>d(P") &e t/ie Veronese 
variety, let 8 = [|J, cmd let r < rs, n - Assume that we are given a point 

P = Pi +P 2 H-h Pd & oy(V). 

Let N be the following linear space: 

N = ker(</>^) o image(^) _L C S' d (C n+1 )*, 

i.e., the symmetric product of the kernel and the complement of the image of </>((. If 
rank</)^ = r, and dimN = ^ ^ ^ — r(n + l), 
then p is a smooth point of a r (vd( P")). 

Proof. The subspace N is the normal space at p of the locus of (r + l)-minors 
of the catalecticant matrix <^; see, for example, [25] Prop. 5.3.3.1]. If N has 
the expected dimension ("J d ) — r(n + 1), then the locus of (r + l)-minors of the 
catalecticant matrix <j>$ is smooth at p and of the expected dimension r(n + 1) — 1. 
The r-secant variety <J r {vd( P")) is contained in that locus, being of the expected 
dimension r(n + 1) — 1 by the Alexander-Hirscliowitz theorem, so it too has to be 
smooth at p. □ 

A natural question concerning the foregoing lemma concerns the maximum value 
of r for which it can be applied. That is, if we pick a sufficiently general smooth 
point p E oy(V), what is the maximum value of r for which Lemma |5.41 can prove 
that p is, indeed, smooth? A lower bound follows immediately from the work of 
Iarrobino and Kanev ,231 Theorem 4.10A]: 

Proposition 5.5. Let V = Ud(P") be the Veronese variety, let 8 = [|J , and let 

n + 8 — 1\ 

j-i y 

Then, Lemma \5.4\ can be applied to all points of an irreducible component of cr r (V) 
minus some Zariski-closed set. 

Proof. The claim follows from [23] and the fact that the conditions on the dimension 
of M and the rank of (jP 5 are valid on dense open sets in the Zariski topology. □ 

In Table [TJ some values of the lower bound in Proposition l5.5l are tabulated along 
with a sharp maximum value of r for which the equations generated by Lemma 15.41 
generate an irreducible component of oy(V). The values of this alleged sharp upper 
bound were computed by taking random points on this variety and verifying Lemma 
15.41 as such, they are only true with high probability. It is clear from the table that 
the lower bound in Proposition 15.51 is not sharp0 

^Recall that the range of applicability of [l8l Proposition 1.32] may be wider by combining it 
with reshapings. 
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Table 1 . The maximum value r for which Lemma 15.41 applies to 
all points in an irreducible component of oy(V) minus some Zariski- 
closed set is displayed as the middle set of columns (X) for each 
degree d = 4, 5,6, 7, 8 of the Veronese embedding V = ^(P”). The 
left set of columns (X) shows the lower bound from Proposition 
1531 for every d. The right set of columns (□) shows the maximum 
value of r for which Kruskal’s criterion is applicable, for every 
d. A * indicates that the value could not be computed within a 
reasonable time. Values displayed in boldface indicate the widest 
range for r for a particular combination of the degree d and size n. 


n 








d 









4 



5 



6 



7 



8 


* 

* 

□ 

* 

X 

□ 

X 

X 

□ 

X 

X 

□ 

X 

X 

□ 

1 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

4 

4 

4 

4 

2 

3 

4 

4 

3 

4 

5 

6 

6 

6 

6 

7 

7 

10 

10 

8 

3 

4 

5 

6 

4 

6 

8 

10 

12 

9 

10 

15 

11 

20 

23 

12 

4 

5 

7 

8 

5 

9 

10 

15 

21 

12 

15 

27 

14 

35 

47 

16 

5 

6 

10 

10 

6 

14 

13 

21 

33 

15 

21 

~k 

18 

56 

87 

20 

6 

7 

12 

12 

7 

19 

15 

28 

50 

18 

28 

~k 

21 

84 

Vc 

24 

7 

8 

16 

14 

8 

25 

18 

36 

72 

21 

36 

~k 

25 

120 

~k 

28 

8 

9 

20 

16 

9 

33 

20 

45 

~k 

24 

45 

~k 

28 

165 

~k 

32 

9 

10 

25 

18 

10 

41 

23 

55 

-k 

27 

55 

~k 

32 

220 

~k 

36 

10 

11 

29 

20 

11 

~k 

25 

66 

-k 

30 

66 

~k 

35 

286 

~k 

40 


Table 2. The maximum value r for which Young flattenings gen¬ 
erate an irreducible component of oy(V) is given in the first row 
(♦). The second row (<0) shows the maximum value of r for which 
the state-of-the-art criterion of PI Proposition f .32] is applicable. 
Values displayed in boldface indicate the widest range for r for a 
particular combination of the degree d and size n. 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

X 

2 

3 

5 

6 

8 

11 

11 

14 

15 

0 

2 

3 

5 

6 

8 

9 

11 

12 

14 


For d = 3, the symmetric flattenings are only sufficient for r = 1 and 2. One 
should employ Young flattenings [55] for extending the range instead. As an illus¬ 
tration of the range that can be covered by such equations, we present in Table 
[2] the maximal value of r for which the equations generated by Young flattenings 
generate an irreducible component of oy(t; 3 (P n )). The range of applicability is also 
compared with the criterion of Domanov and De Lathauwer m Proposition 1.32], 
which to the best of our knowledge provides the state-of-the-art specific identifiabil- 
ity criterion. The border case oq^t^P 6 )) in Tablc[2]is special, because in all of the 
random choices of p £ 011(^3 (P 6 )) that we tested, the Young flattening had rank 
two less than expected. Still, the corresponding minors of the Young flattening did 
cut cr 11 ( 1)3 (P 6 )) scheme-theoretically at p. 
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An implementation in Macaulay2 of the above sufficient condition for smoothness 
based on a symmetric flattening is included in the specif ic-identif iability.m2 
file that is provided with the arXiv version of this paper. In addition, this file 
contains an implementation of a smoothness test based on Young flattenings for 
Veronese embeddings of degree three. 

5.3. An elementary algorithm. For the sake of completeness, we present an 
algorithm that attempts to prove the identifiability of a given Waring decomposition 

p = pi+p 2 ~\ - \-p r £ oy(V), where V = Ud(P”), by checking the sufficient conditions 

in Lemma l5Jl and Lemma PT~il It operates as follows. 

51. Construct a matrix representation of the span M = (T P 1 V, T P2 V,..., T Pr V). 
If rankM < r(n+l), then the algorithm terminates, claiming that it cannot 
prove the identifiability of p. 

52. Construct the symmetric flattening cj) p for 5 = [fj- If rank<^ < r, then 
the algorithm terminates, claiming that it cannot prove identifiability of p. 

53. Compute the matrix N = ker(</>^)oimageIfrankN > ( n ^ d ) —r(n+l), 
then the algorithm terminates, claiming that it cannot prove identifiability 
of p. 

54. Compute a basis of the kernel of M. Denote the number of equations by i. 

55. For every point Pi, i = 1, 2,..., r, perform the following: 

S5a. Construct the Hessians H k for k = 1, 2,... ,£ evaluated at the point 
Pi and stack them into the matrix H. 

S5b. If rankIL < n, then the contact locus is of positive dimension at pi. 
The algorithm halts, claiming that it cannot prove identifiability of p. 

5 6 . The algorithm proclaims that the Waring decomposition p = pi + ■ ■ ■ + p r 
is unique and that p is a smooth point of oy(V). 

It is instructive to investigate the largest value of r for which the above algorithm 
may be expected to prove identifiability of a sufficiently general point p £ oy(V) on 
a generically identifiable Veronese variety V. As the Hessian criterion applies for all 
tensors of subgeneric rank, it follows that the range of applicability is bounded only 
by the smoothness test. That is, the highlighted columns in Table |T| contain the 
relevant values. A popular criterion for testing identifiability that is applicable for 
Veronese embeddings of any degree d > 3 is the so-called Kruskal condition ( 11131 - 
Let pi = af d € V be some specific points. In the symmetric setting, Kruskal’s 
condition states that if 

r < — (dk — d A 1), 

where k is the largest number such that every subset of {a, },; consisting of k vectors 
is linearly independent. For points in general configuration, the maximum value 
for k is thus n + 1. A comparison between the proposed criterion for specific 
identifiability and Kruskal’s criterion is also featured in Table [1] 

5.4. The algorithm at work for a specific example. Consider the following 
10-term Waring decomposition in S 3 C 7 : 

P = J2x 3 i + 

i =0 

(4xo A 3a;i A 2 x 2 A X 3 A 2 x 4 A 3 x 5 A 4xg) 3 A 
(xo A xi A 2 x 2 A 2 x 3 A 3 x 4 A 3x 5 A 4x 6 ) 3 
(xq A 2xi A 3x 2 + 4x 3 A 5 x 4 A 6 x 5 A 7x 6 ) 3 . 
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The identiffability of this example cannot be handled with Domanov and De Lath- 
auwer’s criterion, because it only applies for Waring decompositions with r < 9. 
The sufficient condition presented in this paper, on the other hand, is applicable up 
to 11 terms. Therefore, we can run the algorithm presented above. This example 
may be verified with the specific-identifiability.m2 script that is provided 
with the arXiv version of this paper. In the first step, the 84 x 70 matrix repre¬ 
senting the span of M is constructed. Its rank is 70, as expected. The algorithm 
proceeds with the construction of Young flattening, checking they generate an ir¬ 
reducible component. Then, a basis of the kernel of M is computed, containing 14 
equations. Note that the dimension of N and the codimension of M should always 
be equal in the approach for certifying identifiability that was proposed in this sec¬ 
tion. Then, for each of the points, the 7x7 Hessian matrices are computed and 
stacked into a 7 x 98 matrix H. The rank of H equals 6 for each of the points. 
In addition, the kernel of H T consists of a single vector that must be a multiple 
of the coefficient vector of pp, for example, the vector in the kernel of the stacked 
Hessian H corresponding to the last term in the Waring decomposition of p is a 
multiple of its coefficient vector [ 1234567 ]. For each of the points, this is indeed 
the case. Finally, the algorithm positively concludes that p admits a unique Waring 
decomposition, i.e., p is identifiable. 
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